Vibron-polaron critical localization in a finite size molecular nanowire 
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La Bouloie, 



The small polaron theory is applied to describe the vibron dynamics in an adsorbed nanowire with 
a special emphasis onto finite size effects. It is shown that the finite size of the nanowire discrimi- 
nates between side molecules and core molecules which experience a different dressing mechanism. 
Moreover, the inhomogeneous behavior of the polaron hopping constant is established and it is 
shown that the core hopping constant depends on the lattice size. However, the property of a lattice 
with translational invariance is recovered when the size of the nanowire is greater than a critical 
value. Finally, it is pointed out that these features yield the occurrence of high energy localized 
states which both the nature and the number are summarized in a phase diagram in terms of the 
relevant parameters of the problem (small polaron binding energy, temperature, lattice size). 
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PACS numbers: 



I. INTRODUCTION 



As a result of following Moore's law during the last 40 
years the industry now manufactures integrated cir- 
cuits with feature sizes bellow 180 run. To gain one or 
two orders of magnitude in these feature sizes constitutes 
one of the challenges of the modern technology. In that 
context, among the different ways investigated to con- 
trol the information transfer at the nanoscale, the use of 
adsorbed nanostructures appears promising 0, 0, 0, • 
Indeed, it is now well established that surfaces exhibiting 
self-organized defects are ideal templates for the forma- 
tion of low dimensional nanostructures such as ID wires 
and confined 2D monolayers 0, S 13 • In addition, the 
interest in adsorbed nanodevices is reinforced by the fact 
that local probes, such as scanning tunneling microscope 
(STMV can serve as tools to design the nanostructure 
0, 0, 0, 0, and to excite its electronic or vibra- 
tional degrees of freedom 0, . 

Although the electron transport is expected to play a 
key role in such nanodevices |16L fl7j| . the ability to excite 
the vibration of the admolecules with a STM were ex- 
ploited in recent works to propose an alternative based 
on the use of vibrational excitons |18t fl9t |20| . Indeed, 
in a molecular lattice, the lateral interactions induce a 
coupling between the internal degrees of freedom of the 
molecules and favor the coherent propagation of the in- 
ternal vibrations from one site to another. The dereal- 
ization of the vibrational energy yields the formation of 
vibrational excitons called vibrons. 

In a nanowire, the vibrons do not propagate freely but 
interact with their surrounding formed by the low fre- 
quency modes of the adsorbates and by the phonons and 
the electrons of the substrate. All these excitations are 
responsible for both energy and phase relaxation, which 
the knowledge is essential to understand and exploit the 
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vibronic transport. Note that using the vibrons as a vehi- 
cle for the information transfer requires their confinement 
inside the adsorbate so that metallic surfaces cannot be 
used. Indeed, fast energy relaxation takes place because 
of the resonances of the vibrons with the electron-hole 
pair continuum of the metal which leads to a shorten- 
ing of the vibron lifetime, typically about 1 ps [23|. As 
a result, we restrict our attention onto the use of dielec- 
tric surfaces for which energy relaxation occurs via multi- 
phonon excitations, only, leading to a longer lifetime. For 
instance, the vibrational lifetime for the Si-H vibration 
for the system Si-H is about 1 ns at room temperature 
[22H23| . A remarkable example is given by the stretching 
vibration of CO adsorbed on NaCl for which the lifetime 
was measured to be about 0.01 s 0|- 

In that context, most of the theoretical analysis of the 
vibron dynamics were applied to interpret the infrared 
spectra with a special emphasis on the processes respon- 
sible for line broadening [liliillllliliSllElllllllsl 
|34| . Dephasing mechanisms were thus accounted by as- 
suming a sufficiently weak vibron-phonon coupling and 
by invoking the Markovian limit. Within the standard 
perturbative theory, the effect of the bath was summa- 
rized in a simple parameter : the so-called dephasing 
constant. 

In the present paper, we point out that these as- 
sumptions are questionable to treat the vibron dynam- 
ics in a nanoscale molecular adsorbate. Indeed, it is 
well known that such systems exhibit rather small vibron 
bandwidths, typically about a few wave numbers p2ll33|. 
Therefore, it seems reasonable to assume that the vibron 
hopping constant is about or less than the strength of 
the vibron-phonon coupling so that the standard pertur- 
bative theory cannot be applied. In addition, the time 
spends by a vibron to propagate coherently along one 
nanometer is about a few picoseconds, i.e. the same or- 
der of magnitude than the admitted value for the thermal 
bath correlation time. As a consequence, non Markovian 
effects are expected to play a key role for nanoscale trans- 
port. 
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To overcome these difficulties, we introduce a small po- 
laron theory to describe the vibron transport. In analogy 
with the vibrational dynamics in model proteins , the 
vibron dynamics is s upp osed to be governed by the so- 
called dressing effect HH EE EH E3, Eg] . When 
the vibron-phonon interaction is sufficiently strong, the 
creation of a vibron is accompanied by a lattice distor- 
tion. However, since the vibron bandwidth is smaller 
than the phonon cutoff frequency, the non-adiabatic limit 
is reached 39] . The lattice distortion follows instanta- 
neously the vibron during its propagation so that the vi- 
bron is dressed by a virtual cloud of phonons and forms 
the small polaron. Although the dressing effect is well 
characterized in a lattice with translational invariance, 
very little is known about finite size effects. Therefore, 
in the present paper, the small polaron theory is applied 
to an adsorbed nanowire with a special emphasis onto the 
symmetry breaking induced by the vibron-phonon con- 
finement. Note that the present paper is devoted to the 
characterization of the polaron eigenstates, only. The 
study of the vibron dynamics versus time will be ad- 
dressed in a forthcoming work. 

The paper is organized as follows. In Sec. II, the 
vibron-phonon Hamiltonian used to described the finite 
size nanowire is introduced. In Sec. Ill, the nature of the 
phonon states in a confined nanowire is first summarized. 
Then, a Lang-Firsov |4^| transformation is applied to 
renormalize the vibron-phonon interaction and to reach 
the small polaron point of view. Finally, a mean field pro- 
cedure is performed to obtain the effective Hamiltonian 
for the dressed vibrons. In Sec. IV, a detailed analysis of 
the polaron eigenstates is performed and the results are 
discussed and interpreted in Sec. V. 

II. MODEL AND HAMILTONIAN 

Let us consider a set of N molecules adsorbed on the 
surface of a well-organized substrate. The molecules form 
a finite size ID lattice which the site position is denoted 
n = l,2,...,JV. 

To characterize the vibron dynamics, we assume that 
each molecule n behaves as an internal high frequency 
oscillator described by the standard vibron operators 
and b n . The vibron Hamiltonian is thus expressed as 
(using the convention H = 1) 

N N-l 

H v = £ ^blbn + £ J[bib n+1 + bl +1 b n ] (1) 

n— 1 n— 1 

where luq stands for the internal frequency of each 
molecule and where J denotes the vibron hopping con- 
stant between nearest neighbor admolecules. Note that 
long range lateral interactions may affect the vibron dy- 
namics but these effects are expected to be rather weak 
in a ID lattice. 

To mimic the influence of the phonon bath in a sim- 
ple way, we assume that the internal vibration of each 



molecule interacts preferentially with the longitudinal 
phonons of the nanowire. These phonons refer to the col- 
lective dynamics of the center of mass of the admolecules 
around their equilibrium positions. Within the harmonic 
approximation, the phonon Hamiltonian is thus written 
as 

N v 2 N ~ l 1 

h p = £ ^ + £ 2 w[Un+l " Un)2 (2) 

n— 1 n—1 

where u n and p n denote the displacement and the corre- 
sponding momentum of the n th admolecule, respectively. 
M stands for the mass of each admolecule which interact 
with its nearest neighbor molecules via the lateral force 
constant W. 

By neglecting the population relaxation, the phonon 
bath leads to stochastic fluctuations of the dynami- 
cal parameters which characterize the vibron dynamics. 
In this work, we use the potential deformation model 
EE EE EE E3 which supposes that the bath induces 
a random modulation of the internal frequency of each 
molecule. Within this model, the vibron-phonon cou- 
pling Hamiltonian is expressed as 

AT-l 

AH vp = 22 x(u n +i - u„_i)&t b„ 

n=2 

+ x(«2 - ui)b\bi + x(v>n - u N ^i)b' N b N (3) 

where x IS the strength of the vibron-phonon coupling. 
Note that the previous model cannot be applied to 
molecules on metals. Indeed, Persson and co-workers (see 
for instance [2H l28l [2^) have shown that the high fre- 
quency stretching mode is preferentially coupled to either 
the frustrated translational motion or the frustrated rota- 
tional motion of the admolecule. Since the admolecules 
are adsorbed onto high symmetry adsorption sites, the 
vibron-phonon coupling depends in a quadratic way on 
the phonon coordinates. However, as pointed out in the 
introduction, we restrict our attention to admolecules 
weakly bounded to the surface for which the preferen- 
tial pathway for dephasing involves the collective modes 
of the monolayer. This situation occurs, for instance, for 
the CO/NaCl system 0. 

Finally, the vibron-phonon dynamics in the nanowire 
is governed by the full Hamiltonian H = H v +H p +AH vp . 
The following section is devoted to its simplification and 
to the introduction of an effective Hamiltonian for the 
vibrons, only. 

III. THEORETICAL BACKGROUND 

As mentioned in the introduction, the vibron dynam- 
ics in a molecular lattice with a strong vibron-phonon 
interaction is essentially governed by the dressing effect. 
The characterization of this effect requires the knowledge 
of the phonon eigenstates and of the Lang-Firsov trans- 
formation [44j which yields the polaron point of view. 
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However, although these two features are well known in 
an infinite lattice with translational invariance, finite size 
induces a symmetry breaking responsible for a modifica- 
tion of both the phonon structure and the dressing mech- 
anism. The present section is thus devoted to the pre- 
sentation of these modifications. 



A. phonon states 

In a finite size nanowire, the phonon states do not cor- 
respond to Bloch waves with well defined wave vectors. 
Indeed, the sides of the lattice lead to a reflexion of the 
acoustic waves so that the true eigenstates appear as a 
superimposition of incident and reflected plane waves. 
More precisely, the nanowire, which is responsible for the 
confinement of the phonons, behaves as a resonant cav- 
ity. It selects particular values of the wave vector due to 
the free boundary conditions so that a stationary regime 
takes place. 

Therefore, it is straightforward to show that the 
phonon eigenstates correspond to N normal modes with 
quantized wave vectors q p = pn/N, p = 0,l,..,N—l and 
frequency f2 p = Q c sin(p7r/2iV), where f2 c = y/4W/M 
denotes the cutoff frequency of the corresponding infi- 
nite nanowire. Within this normal mode representation, 
the phonon Hamiltonian H p (Eq.JSJ) is thus rewritten as 



N—\ 

H P =J2 n P (<4> a P + V 2 ) 

p=0 



(4) 



The vibron-phonon coupling Hamiltonian AH vp (Eq.(|3J)) 
is expressed as 

N JV-l 

AH vp = ^ Xp,n( a p + a p) b l b n (5) 

n— 1 p— 1 

where Xp.n characterizes the strength of the coupling be- 
tween the vibrons and the p th normal mode. From 
Eqs.JHJl, it is straightforward to show that 



Xp,: 



2_ smjpn/N) _ gWg r ( „_ /9)} (g) 
N Vl sin (W2A0 I N 



where Ao = xQ^MW)^ 1 / 4 (ft as be reintroduced to 
avoid confusion) . Note that the mode with zero wave vec- 
tor describes the free translation of the nanowire. How- 
ever, since this mode disappears in an adsorbed nanowire 
due to its interaction with the substrate, it will be disre- 
garded in the following of the text. 



B. Small polaron theory and effective Hamiltonian 

To partially remove the vibron-phonon coupling 
Hamiltonian, a Lang-Firsov transformation is applied 



|44j . Indeed, since the vibron-phonon dynamics is domi- 
nated by the so-called dressing effect, we consider a "full 
dressing" and introduce the following unitary transfor- 
mation 



JV N-l 

U = exp(£ - a vP n b n ) (7) 

n—1 p— 1 " 

By using Eq.JJJ, the transformed Hamiltonian H = 
U HU^ is written as 



JV 



n=l 

JV-1 

+ Y A^lfin+iblfin+i + h.o.} - Esb'l^bibn^bn 

n=l 
N-l 

+ Y n P a l a P ( 8 ) 



where h.o. denotes the hermitian operator and where 
e„ = Eb(1 — l/2((5 ri: i + <5 71j iv)) is expressed in terms of the 
small polaron binding energy Eb = 2Aq/H c . In Eq.©, 
G„ stands for the dressing operator defined as 



6„ = exp(- 



N-l 

Xn, P t 



dp)) 



(9) 



0=1 



p 



In this dressed vibron point of view (Eq.JSJ)), the 
vibron-phonon coupling remains through the modulation 
of the lateral terms by the dressing operators. Although 
these operators depend on the phonon coordinates in a 
highly nonlinear way, the vibron-phonon interaction has 
been strongly reduced within this transformation. As a 
result, we can take advantage of such a reduction to per- 
form a mean field procedure 38, 39] and to express the 
full Hamiltonian H as the sum of three separated contri- 
butions as 



H = H. 



eff 



H V + AH 



(10) 



where H e ff = ((H — H p )) denotes the effective Hamilto- 
nian of the dressed vibrons and where AH = H — H p — 
((H — H p )) stands for the remaining part of the vibron- 
phonon interaction. The symbol (. . .} represents a ther- 
mal average over the phonon degrees of freedom which 
are assumed to be in thermal equilibrium at temperature 
T. 

As a result, the effective dressed vibron Hamiltonian 
is written as 

N 
ri=l 

N-l 

7 (,+ 1 

+ h.o.] — £ , B 6^ +1 bJ l 6„ + i6„ (11) 
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where J n = Jexp(— S(T, n)) and where S(T,n) is the 
coupling constant defined as (/«b denotes the Boltzmann 
constant) 

1 JV_1 - O 
S(T, n) = ±Y, ( *"'" n Xn+1 ' P ) 2 c ° th (-^) (12) 
p—i p ^ 

The Hamiltonian H e ff fEa.l|lip) describes the dynam- 
ics of vibrons dressed by a virtual cloud of phonons, i.e. 
small polarons. It accounts for a renormalization of the 
main part of the vibron-phonon coupling within the non- 
adiabatic limit. The remaining vibron-phonon coupling 
is thus assumed to be small in order to be treated using 
perturbative theory. Such a contribution, disregarded in 
the present work, will be addressed in a forthcoming pa- 
per. 

Ea. Hll|) clearly shows the interplay between the usual 
dressing effect and the finite size of the nanowire. In- 
deed, within our procedure, the main features induced 
by the dressing effect in an infinite lattice are recovered. 
More precisely, the dressing modifies the harmonic part 
of the Hamiltonian and yields a redshift of each internal 
frequency. Moreover, it reduces the strength of the vi- 
bron hopping constant so that the effective mass of the 
polarons increases. In addition, the dressing favors the 
occurrence of coupling terms which break the vibron in- 
dependence and which directly affect the multi-vibron 
dynamics 0> C3- These terms will be disregarded in 
the present work devoted to the single- vibron dynamics. 

However, the finite size of the lattice discriminates be- 
tween the core molecules, i.e. those which are located 
at the sites n — 2, ...,7V — 1, and the side molecules 
n = 1 and n = N. Indeed, Eq.Q shows that core 
and side molecules do not experience the same inter- 
action with the phonon bath which results in a differ- 
ent dressing effect. The redshift of the frequency of the 
core molecules, equal to the small polaron binding en- 
ergy Eb, corresponds to the shift occurring in an infinite 
lattice. In a marked contrast, the side molecules experi- 
ence a frequency shift two times smaller. In other words, 
for the polaron dynamics, the finite size nanowire ex- 
hibits two defects located at the side sites n = 1 and 
n = N . In addition, the modification of the nature of 
the phonons induced by their confinement favors an in- 
homogeneity in the polaron hopping constant. Indeed, 
in a marked contrast with the translational invariant lat- 
tice, the coupling constant S(T,n) depends on the site 
position. As a result, the effective hopping constant 
J n = Jexp(— S(T, n)) is not uniform along the nanowire 
and enhances the symmetry breaking induced by the con- 
finement. 

At this step, the restriction of the effective Hamilto- 
nian H e ff (Eq.l^J)) to the single- vibron subspace can 
be easily diagonalized by using standard numerical pro- 
cedures. As a result, the behavior of the single-vibron 
eigenstates can be characterized depending on the values 
of the parameters involved in the model as illustrated in 
the following section. 



IV. NUMERICAL RESULTS 

In this section, the previous formalism is applied to 
characterize the small polaron eigenstates in a finite size 
nanowire. The undressed hopping constant is fixed to the 
typical value J — 4.0 cm -1 whereas the phonon cutoff 
frequency is set to J7 C = 100 cm^ 1 . Note that the fre- 
quencies are centered around the vibrational frequency 
wo set to zero. The temperature T, the small polaron 
binding energy Eb and the lattice size N are parameters 
allowed to vary. 

In Fig. 1, the dependence of the polaron hopping con- 
stant J n on the position of the sites is drawn for T = 5 
K (Fig. la) and T = 50 K (Fig. lb). The num- 
ber of molecules is fixed to N — 20 and three typical 
values for Eb are considered. Whatever the temper- 
ature, the hopping constant decreases as Eb increases 
as a result of the well-known dressing effect. At low 
temperature (Fig. la), the hopping constant displays 
a rather important inhomogeneity which increases as Eb 
increases. Although J n is almost uniform in the core 
of the nanowire, it increases close to the lattice sides 
where a reduced dressing effect takes place. For weak 
Eb values, i.e. Eb = 10 cm -1 , only the hopping con- 
stant involving side molecules differs significantly from 
the hopping constant in the core of the nanowire. For 
instance, the ratio Ji/J — 0.95 is slightly greater than 
the core ratio equal to Jn/2/J = 0.92. By contrast, as 
when increasing Eb, both the strength and the range 
of the inhomogeneity increase. When Eb = 100 cm" 1 , 
the ratio Ji/J = 0.58 is 1.3 times greater than the core 
ratio equal to — 0.43 and the hopping constant 

J2/J = 0.45 slightly differs form the core value. At high 
temperature (Fig. lb), the inhomogeneous effects are 
drastically reduced. Only the hopping constant involv- 
ing side molecules differs form its core value whatever 
Eb- However, this difference represents less than 10% of 
the core hopping constant. Note that as when increasing 
again the temperature, the inhomogeneous effects almost 
disappear. 

To characterize the inhomogeneous nature of the hop- 
ping constant, Fig. 2 represents the behavior of the dif- 
ference between the side and the core hopping constants 
with respect to E B for N = 30 and T = 5 K (full cir- 
cles), 50 K (open circles) and 100 K (full triangles). The 
illustrated features coroborate the previous results since 
it is shown that the difference between the side and the 
core hopping constants increases as Eb increases and as 
the temperature decreases. 

Finally, the behavior of the core hopping constant with 
respect to the size N of the nanowire is illustrated in Fig. 
3 for T = 50 K and for E B = 10 cm" 1 (full circles), 
50 cm -1 (open circles) and 100 cm -1 (full triangles). 
As previously, the decrease of the core hopping constant 
with the small polaron binding energy is a signature of 
the dressing effect. However, Fig. 3 clearly shows that 
the dressing mechanism is size dependent and the core 
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hopping constant exhibits two regimes with respect to 
N. For small sizes, typically when N < 20, the core 
hopping constant decreases rapidly as N increases. Note 
that it is equal to the undressed value J when N = 2. By 
contrast, when the lattice size is sufficiently important, 
the core hopping constant converges to a constant value 
equal to the polaron hopping constant in an infinite lat- 
tice. In other words, the dressing effect decreases as the 
size of the nanowire decreases. 

In Figs. 4, 5 and 6, the polaron cigenstates of a 
nanowire with TV = 15 molecules are characterized for 
three different Eb values. The temperature is fixed to 
T = 100 K so that the inhomogeneity of the polaron 
hopping constant affects the sides of the nanowire, only. 
We thus consider that the hopping constant is uniform 
in the core of the lattice and equal to its value Jn/2 & t 
the center of the nanowire. Such a behavior allows us 
to introduce the polaron band which the eigenenergies ui 
range between ujq — Eb — 2Jn/2 an d — Eb + 2J N / 2 - 
This band contains the single polaron states of an in- 
finite nanowire, i.e. Bloch waves, characterized by the 
finite size hopping constant Jn/2- Therefore, to charac- 
terize this band for the three different Eb values, the 
corresponding densities of states are drawn (full line) in 
Figs. 4a, 5a and 6a. 

When Eb = 3 cm -1 (Fig. 4a), the core hopping con- 
stant is equal to Jn/2 — 3.71 cm -1 and the polaron band 
extends from —10.42 cm -1 to 4.42 cm -1 (around uo ). 
The finite size nanowire exhibits N = 15 eigenstates 
which the energies are located inside the band. To il- 
lustrate the nature of these states, the wave function of 
the two states with the two higher energies are drawn in 
Fig. 4b. The figure clearly shows that these states cor- 
respond to the first two stationary states of a finite size 
lattice. In fact, in that case, our calculations establish 
that the N eigenstates are N stationary states. 

When Eb = 15 cm -1 (Fig. 5a), the polaron band- 
width is reduced (Jn/2 = 2.74 cm -1 ) and extends from 
—20.48 cm -1 to —9.52 cm -1 . The nanowire supports 
N — 2 = 13 states which the energies are located inside 
the band and exhibits two states outside the polaron 
band. The corresponding energies, almost degenerated 
and equal to —6.49 cm -1 , are greater than the top of 
the band. As shown in Fig. 5b, these two states corre- 
spond to two localized states which the wave function is 
important close to the sides of the nanowire. Note that 
one state appears symmetrically localized onto each side 
whereas the other one exhibits an anti-symmetric local- 
ization. 

As shown in Fig. 6a, an intermediate situation occurs 
when Eb = 7.2 cm -1 . In that case, the polaron band 
extends from —13.86 cm -1 to —0.54 cm -1 which cor- 
responds to a core hopping constant Jn/2 — 3-33 cm -1 . 
Our calculations clearly show that a single cigenstate has 
an energy equal to —0.48 cm -1 greater than the top of 
the band. A shown in Fig. 6b, this state corresponds to 
a state symmetrically localized onto each side. However, 
the localization length is larger than or about to the sys- 



tem size so that the wave function exhibits a significant 
value over the entire lattice. 

The previous results illustrate the fact that the con- 
fined nanowire exhibits localized states depending on the 
values of the small polaron binding energy. In fact, our 
calculations have established that the finite size lattice 
supports at least zero, one or two localized states. To 
characterize this feature, the behavior of the energy of the 
localized states vs Eb is displayed in Fig. 7 for different 
N values and for T = 100 K. Note that the origin of the 
energy is taken at the top of the polaron band. The figure 
clearly shows that the number of localized states depends 
on the value of both Eb and N. When N = 10, the first 
localized state occurs when Eb = 6.79 cm -1 (full cir- 
cles) whereas the second localized state takes place when 
Eb = 8.06 cm -1 (full triangle). The energy of each local- 
ized state scales similarly with respect to the small po- 
laron binding energy. More precisely, the energy exhibits 
a power law dependence close to the top of the band and 
reaches a linear regime for sufficiently strong Eb values. 
In this linear regime, the energies of the two localized 
states are almost degenerated. When N = 15, the first 
localized state occurs when Eb = 6.71 cm -1 (open cir- 
cle), i.e. almost the same Eb value responsible for the 
occurrence of the first localized state when N = 10. By 
contrast, the increase of the size of the nanowire reduces 
the critical Eb value for the occurrence of the second lo- 
calized states. Indeed, when N = 15, this latter state 
takes place when Eb = 7.51 cm -1 (open triangles). 

In the same way, Fig. 8 shows the behavior of the 
energy of the localized states versus the temperature for 
different N values and for Eb = 7 cm -1 . The figure 
clearly shows that the nanowire supports zero localized 
states at low temperature. However, as when increas- 
ing the temperature, two transitions appear successively 
leading to the occurrence of two localized states. When 
N = 10, the first localized state occurs at T = 50 K (full 
circles) whereas the second localized state takes place at 
T = 146 K (full triangle). The energy of each localized 
state scales similarly with respect to the temperature. It 
first exhibits a power law dependence close to the top of 
the polaron band and reaches a linear regime at a suffi- 
ciently high temperature. In the linear regime, the en- 
ergies of the two localized states are degenerated. When 
N = 15, the first localized state occurs at T = 45 cm -1 
(open circles) whereas the second localized state is cre- 
ated at T = 106 K (open triangles). Note that as when 
increasing the value of the small polaron energy, the crit- 
ical temperatures at which the localization takes place 
are reduced. 



V. DISCUSSION 

To interpret and discuss the previous numerical results, 
let us first focus our attention onto the nature of the 
polaron hopping constant. Indeed, it has been shown (see 
Figs. 1-3) that the finite size of the nanowire favors an 
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inhomogeneous behavior of the hopping constant which 
appears site dependent. However, such a feature occurs 
essentially at low temperature and strong small polaron 
binding energy and the hopping constant becomes almost 
site independent at high temperature. 

Indeed, in this latter situation, it is straightforward 
to show that the coupling constant S(T,n) Ea. i|12|) is 
expressed as 



S(T,n) = — ^— (1--) 



(13) 



As a consequence, the effective hopping constant within 
the high temperature limit is written as 



(14) 



where denotes the polaron hopping constant in an 
infinite nanowire and where TV* cx EsksT /Vl 2 c . Eq. ifH|) 
shows that the high temperature polaron hopping con- 
stant is translationaUy invariant but depends on the size 
of the nanowire. More precisely, it increases as the lat- 
tice size decreases to reach its undressed value equal to 
J when the nanowire reduces to a dimer with N = 2. In 
that case, the phonon bath induces the same frequency 
modulation for the two molecules n = 1 and n = 2 (see 
Eq.© when N — 2). Therefore, the resonance between 
the two internal modes is never broken by the phonon 
motion so that nothing prevents the derealization of the 
vibrational energy. The vibrational energy flow occurs 
according to the rate J which characterizes undressed 
molecules. By contrast, as when increasing the lattice 
size, the polaron hopping reach the limiting value Joo of 
the effective hopping constant in an infinite lattice. 

In other words, the dressing effect is size dependent. 
This feature originates in the modification of the phonon 
states induced by their confinement in the finite size 
nanowire. Indeed, this confinement is responsible for the 
occurrence of a finite number of stationary states char- 
acterized by a discrete energy spectrum. As decreasing 
the lattice size, the number of phonon modes decreases so 
that the vibron-phonon interaction is reduced when com- 
pared to that in a lattice with the translational invari- 
ance. However, our results clearly show that the prop- 
erty of the polaron band is recovered when the size of the 
nanowire is greater than the typical value N* (see Fig. 
3). 

At low temperature, no analytical expression of the 
coupling constant can be obtained. However, we can take 
advantage of the fact that the dependence of the hopping 
constant on the site position is important in the vicinity 
of the nanowire sides, only (see Fig.l). Therefore, to 
described such a behavior, we consider a semi-infinite 
nanowire and perform the calculation of Ea. l|12l) within 
the limit N — > oo. As a result, at zero temperature, the 
coupling constant near the side n = 1 is expressed as 



S(0,n) 



16E B 9 - 104n 2 + 128n 4 
3tt^ c 9 - 160n 2 + 256?i 4 



Ea.pSjl allows us to reproduce the behavior of the po- 
laron hopping constant in the vicinity of the side n = 1. 
It shows that the range of the inhomogeneity of the hop- 
ping constant depends on the ratio Eb/Q c - In other 
words when Eb is smaller than the phonon cutoff fre- 
quency the inhomogeneous effects disappear as shown in 
Fig. 1. 

The previous features allow us to characterize the 
single-polaron dynamics in a simple way in order to un- 
derstand the mechanisms involved in the occurrence of 
localized states. Indeed, the numerical results clearly 
show that the inhomogeneity of the hopping constant 
affects the sides of the nanowire, only. Therefore, we 
suppose that the hopping constant is uniform in the core 
of the lattice and equal to its value Jn/2 a t the center 
of the nanowire. By contrast, we assume that hopping 
processes involving side molecules is singular and charac- 
terized by the constant J\ — Jn~i so that the hopping 
constant is approximated as 



J NJ2 if n = 2,..., N -2 
Ji if n = 1 and n = N — 1 



(16) 



As a result, by inserting Ea. (|16fl into the expression of 
the effective polaron Hamiltonian H e ff Eci. (|ll|) yields 
the Schrodinger equation H e ff \ i>) = w | V*) f° r single- 
polaron states as 

Xipi = ocipi + (1 + /3)V>2 

AV>3 = ^2 + 1p4 



\ip N = aip N + (1 + 0)ipN-i 



(17) 



(15) 



where ip n denotes the polaron wave function onto the n 
th site and where A = (lu—luq+Eb)/Jn/2, ct = Eb/2Jn/ 2 

and ft = (Ji - J n /2)/Jn/2- 

From a physical point of view, the wave function ip n is a 
superimposition of Bloch waves ip n = A + e lkn + A~e~ lkn 
that propagate in both directions along the nanowire. 
However, because of the symmetry breaking induced by 
the confinement, the wave vector k is quantized accord- 
ing to N values which can be real or complex, depending 
on the various parameters. A real solution for k char- 
acterizes a stationary wave with a frequency inside the 
polaron band A = 2cos(/c), i.e. ujq — Eb — 2J^/2 < u < 
luq — Eb + 2J/v/2- A complex value of k corresponds to 
a localized state which the wave function is strongly lo- 
calized on the side molecules and which the energy lies 
below or above the polaron band. A localized state is 
characterized by its localization length £, i.e. the inverse 
of the imaginary part of k, which defines the spatial ex- 
tension of its amplitude. 

As detailed in Refs. 0,112, the occurrence of local- 
ized states corresponds to a critical transition between a 
stationary regime and a localized regime as the relevant 
parameters of the problem vary. Therefore, as shown in 
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Appendix A, the Renormalization Group (RG) theory 
can be used in order to understand the transition and to 
predict the occurrence of localized states. 

The conditions for the occurrence of high energy local- 
ized states is thus expressed as 

a + (1 + (3) 2 - 2 = 

a +(l + /5) 2 -2-^- T (l + /3) 2 =0 (18) 

whereas the condition for the occurrence of low energy 
localized states is written as 

a - (1 + (3) 2 + 2 = 

a -(l + /5) 2 + 2+^- T (l + /3) 2 =0 (19) 

As shown in Fig. 9, Eas. ((T%|> and (|T§1) define the crit- 
ical curves in the space of the parameters (a, [3) which 
separate phases exhibiting zero, one, two, three and four 
localized states. More precisely, as indicated in the fig- 
ure, the full lines characterize the occurrence of high en- 
ergy localized states located above the polaron band. The 
lines discriminate between three regimes connected to the 
presence of zero (OHE), one (1HE) or two (2HE) high 
energy localized states. By contrast, the dashed lines 
refer to the occurrence of low energy localized states, lo- 
cated below the polaron band, and discriminate between 
three regimes connected to the presence of zero (OLE), 
one (1LE) or two (2LE) low energy localized states. 

To explain the behavior of the phase diagram, let us 
consider the situation addressed in our model and which 
corresponds to a > and (3 > 0. The generalization 
of the entire parameter space is straightforward. When 
(3 = 0, i.e. when the hopping constant is uniform in the 
nanowire, the lattice exhibits two defects corresponding 
to the side sites which the energy is blueshifted when 
compared to that of the core sites. Therefore, the phase 
diagram shows that two critical values for a discrimi- 
nate between phases with zero, one and two high energy 
localized states. By contrast, when a = 0, the lattice ex- 
hibits two defects corresponding to a blueshift of the hop- 
ping constant which connects each side site to its nearest 
neighbor site. The phase diagram shows that two criti- 
cal values for (3 discriminate between phases with zero, 
two, and four localized states. In the two localized states 
regime, both a low energy and a high energy localized 
state occur whereas in the four localized states regime, 
the system exhibits two low energy and two high energy 
localized states. Finally, for both non vanishing a and (3 
values, intermediate situations take place. 

As detailed in Refs. [&]], 15^, such a behavior can be 
interpreted in terms of a semi-infinite lattice with defects 
at its side. For such a system, one or two localized states 
can occur depending on the values of the parameters. For 
instance, when (3 = 0, one high energy localized state ap- 
pears when a > 2. In the same way, when a = 0, the lat- 
tice exhibits two localized states when (f3 + l) 2 > 2. One 



state is located below the polaron band (low energy local- 
ized state) whereas the second one appears above the po- 
laron band (high energy localized state). The finite size 
nanowire can thus be viewed as the superimposition of 
two semi-infinite lattice with defects at their finite sides. 
Thus, when the required relations are full filled, one (or 
two) localized state occurs on each side of the nanowire. 
However, states localized onto two different sides interfere 
inside the chain core. The strength of the coupling de- 
pends on the overlap of their respective amplitudes which 
leads to hybridization. For a strong overlap, the coupling 
can be sufficient to strongly split the energies and push 
back half of the states into the polaron band which thus 
loses their localized specificity. Therefore, the nanowire 
exhibits one (or two) localized state. By contrast, when 
the overlap is weaker, the energies of the states remain 
outside the polaron band and the nanowire supports two 
(or four) localized states. Since the occurrence of this 
latter regime depends on the overlap of the original lo- 
calized states, the corresponding critical values depend 
on the system size (See Eas. (ll8fl and flU ). For instance, 
when (3 = 0, it is straightforward to show that the semi- 
infinite lattice exhibits a high energy localized state when 
a > 2. Therefore, in the finite size nanowire, two original 
localized states occur on each side when a > 2. How- 
ever, these two states interact due to their overlap inside 
the core of the lattice and their hybridization takes place 
which results in the formation of two states which ap- 
pears as superimpositions of the two original localized 
states. As a result, when 2 < a < 2 + 2/ (N - 1), the hy- 
bridization is strong enough so that the low energy state 
is pushed back into the polaron band. The lattice thus 
supports a single localized state, only. By contrast, when 
a > 2 + 2/(iV — 1), the hybridization is weaker and the 
two states remain above the polaron band. 

Although the phase diagram displayed in Fig. 9 al- 
lows for a complete understanding of the localization 
transition, it involves parameters which are not indepen- 
dent. For instance, both a and [3 depend on Eb and 
on the temperature. However, as illustrated in Fig. 10 
for TV = 10 (full line) and N = 20 (dashed line), we 
can define a phase diagram in the (Eb, T) parameter 
space. In the parameter range considered here, the dia- 
gram exhibits two critical curves which discriminate be- 
tween phases with zero, one and two localized states. The 
figures clearly shows that both the small polaron binding 
energy and the temperature favor the occurrence of local- 
ized states. However, these two parameters do not con- 
tribute in the same way to the localization mechanism. 
The influence of the small polaron binding energy in the 
localization processes is threefold. First, as when increas- 
ing Eb, the energy difference between the side and core 
sites increases. Then, strong Eb values yield a strong 
dressing effect so that the core hopping constant is dras- 
tically reduced. Finally, as when increasing Eb, the inho- 
mogeneous behavior of the polaron hopping constant is 
enhanced (see Fig. 2) . All these features act in the same 
direction and favor the localization. By contrast, the in- 
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fluence of the temperature on the localization processes 
originates essentially in the dressing effect. Indeed, as 
when increasing the temperature the core hopping con- 
stant is reduced and the localization is enhanced. Note 
that the temperature favors a uniform behavior of the 
hopping constant and prevents the localization due to 
the (3 parameter. 

In addition, the phase diagram illustrates the influence 
of the size of the nanowire. It shows that the lower criti- 
cal curve depends slightly on N whereas the upper criti- 
cal curves exhibits a more important dependence. These 
two features have two distinct origins. Indeed, the size 
dependence of the lower critical curve originates in the 
dependence of the core hopping constant with respect 
to the lattice size. By contrast, the size dependence of 
the upper critical curve is essentially due to the overlap 
mechanism involved in the superimposition of the origi- 
nal states localized onto different sides of the nanowire. 
Therefore, Fig. 10 shows that this latter dependence de- 
creases as the size of the lattice is increased, in agreement 
with Eas. (|18|l and l|19(l . so that the two critical curves 
tend to collapse. 

Finally, in the space of the parameters (Eb, T), Fig. 
10 clearly shows that the finite size nanowire can exhibit 
a maximum of two localized states. This feature seems 
to be in contrast with the results displayed in the phase 
diagram connected to the general (a, (3) space where it 
is shown that a maximum of four localized states can be 
reached (see Fig. 9). In fact, since the parameters a 
and P are expressed in terms of Eb and T, the physical 
space {Eb, T) corresponds to a particular projection in 
the general space (a,j3). More precisely, the range for the 
variations of the physical parameters Eb and T generates 
the restriction a > 0, > and a >> (3. In that 
case, Fig. 9 clearly shows that the nanowire supports a 
maximum of two high energy localized states. 

To illustrate these features, let us estimate the order of 
magnitude of the key parameters for the CO/NaCl sys- 
tem from the knowledge of the vibron dynamics in the 
monolayer [3i| . The CO molecules form a (2x1) struc- 
ture with a rectangular unit cell containing two molecules 
mutually perpendicular. Along the longer size of the cell, 
the vibron bandwidth is about 3 cm -1 whereas it reaches 
12 cm -1 along the perpendicular direction. The frequen- 
cies of the external modes of the monolayer range from 
36 cm -1 for the translational motion parallel to the sur- 
face to 140 cm -1 for the collective librational motion. 
The comparison with experimental spectra has shown 
that the phase relaxation is due to the coupling between 
the vibrons and the collective librons of the monolayer. 
Moreover, the dephasing constant exhibits a linear de- 
pendence vs the temperature in agreement with the de- 
formation potential approximation used in the present 
work. Although the small polaron binding energy is un- 
known for this system, we expect a rather weak E# value, 
less than or about to 1 cm -1 , in agreement with the weak 
value of the dephasing constant equal to 0.2 cm -1 at 55 
K. As a consequence, we expect that the weak Eb value 



prevents the formation of any localized state. 

Note that the present model can also be applied to 
characterize the vibron-phonon dynamics in a finite size 
a-helix protein (see for instance Ref. |35|). In that 
case, the high frequency modes correspond to either the 
C=0 stretching vibrations (amide-I) or the N-H stretch- 
ing modes of the peptide groups of the protein. The 
phonons characterize the collective dynamics of the ex- 
ternal motions of the peptide groups. The parameters 
commonly used to described the system are J=7.8 cm" 1 , 
tt c = 87 - 137 cm" 1 , E B =15 cm" 1 for the C=0 vibrons 
(see for instance Refs. |42l l43| p and J=5 cm" 1 , f2 c = 100 
cm" 1 , F,b= 84 cm" 1 for the N-H vibrons |53|. In these 
cases, we expect that the strong Eb value allows for the 
formation of one or two localized states at biological tem- 
perature. Note that the adsorption of a-helices or similar 
molecules could be considered. 



VI. CONCLUSION 

In the present paper, the small polaron theory was 
applied to describe the vibron dynamics in an adsorbed 
nanowire with a special emphasis onto finite size effects. 
It has been shown that side molecules and core molecules 
do not experience the same interaction with the phonon 
bath which results in a different dressing effect. The 
redshift of the frequency of the core molecules is equal 
to the shift occurring in an infinite lattice whereas side 
molecules experience a frequency shift two times smaller. 
In addition, the inhomogeneous behavior of the polaron 
hopping constant has been established. The hopping con- 
stant between a side molecule and its nearest neighbor 
molecule is greater than the hopping constant which con- 
nects core molecules. Moreover, it has been shown that 
the core hopping constant depends on the size of the 
nanowire which indicates that the dressing mechanism is 
size dependent. However, we have shown that the prop- 
erty of a lattice with translational invariance is recovered 
when the size of the nanowire is greater than a typical 
value N* cx EsksT These features have been at- 
tributed to the modification of the phonons induced by 
their confinement which is responsible for the occurrence 
of a finite number of stationary states. As decreasing 
the lattice size, the number of phonon modes decreases 
so that the vibron-phonon interaction is reduced when 
compared to that in a lattice with translational invari- 
ance which results in a softening of the dressing effect. 

The modification of the dressing mechanism was shown 
to be responsible for the occurrence of localized states. 
The nature and the number of localized states versus the 
relevant parameters of the problem have been summa- 
rized in a phase diagram. Although the lattice is able 
to supports a maximum of four localized states, it has 
been shown the for physical values of the parameters, 
the nanowire can exhibit zero, one or two high energy 
localized states. 
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APPENDIX A: DECIMATION OF THE 
SCHRODINGER EQUATION 

To predict the occurrence of localized states, a deci- 
mation 0, Unl E3 01 the Schrodinger equation Eq. 
(|17|l is performed. It consists in eliminating one site over 
two from the initial lattice to arrive at a scaled lattice 
that has twice the initial lattice spacing. In Eqs. 1)170. 
we eliminate the amplitudes of the even sites by substi- 
tuting their expressions in the Schrodinger equations for 
the odd sites and obtain a new set of equations for the 
odd sites only, as 

(A 2 - 2)tpi = (aA + (l + /5) 2 -2)^ 1 + (l + /5)^ 3 (Al) 

(A 2 -2)V> 3 = (l + /3)Vi+V>s 

(A 2 -2)Vs = ^3+^7 

(A 2 - 2)ip N = (aX + (1 + I3) 2 - 2)4> N + (1 + 0)r/> N _ 2 

Eq.(Al) characterizes the Schrodinger equation of a 
rescaled lattice with new parameters defining the RG 
transformation R as 

A (1) = A 2 -2 (A2) 
= aA + (l + /3) 2 -2 

When the decimation is applied recursively p times, with 
N = 2 P+1 + 1, we obtain an ultimate scaled lattice formed 



by three sites, only, for which the eigenenergies are ex- 
pressed as 

A (p) = (A3) 

/ x a (p) RcM 
A (p) = ^-±^/^-) 2 + 2(l+/3(rt)2 

At this step, a high energy localized state occurs when 
its energy is at least equal to the maximum of the po- 
laron band. This condition is obtained when A = A c = 2. 
At the critical point A = A c , the scaled values of the pa- 
rameter A (Eq. (A2)) satisfy A^ = ... = A^ = A c . This 
parameter A becomes scale invariant and is a fixed point 
A c of the RG transformation R that verifies A c = R(X C )- 
Indeed, no change in the critical parameters accompanies 
the length scaling so that the scaled lattice remains at a 
critical point. 

From Eqs .(A3), the scaled values of the dynamical 
parameters at the critical point are expressed as 

a« = 2-(l + /3) 2 -^^(2- a -(l + / 3) 2 XA4) 

/3 (p) = [3 

Combining Eqs. (A3) and (A4) for A (p) = 2 leads to the 
conditions for the occurrence of high energy localized 
states Eas. (|18|) . The condition (|19|l for the occurrence 
of low energy localized state can be obtained in a similar 
way. 
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FIG. 1: Polaron hopping constant J n versus the site position for T = 5 K (a) and T = 50 K (b) and for different values of the 
small polaron binding energy Eb- The number of molecules belonging to the nanowire is fixed to N = 20. In all the figures, 
the undressed hopping constant is fixed to the typical value J = 4.0 cm" 1 and the phonon cutoff frequency is set to fi c = 100 
cm -1 . 
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FIG. 2: Behavior of the difference between the side and the core hopping constants vs the small polaron binding energy Eb 
for N = 30 and T = 5 K (full circles), 50 K (open circles) and 100 K (full triangles). 
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FIG. 3: Behavior of the core hopping constant vs the size N of the nanowire size for T — 50 K and for Eb = 10 cm 1 (full 
circles), 50 cm -1 (open circles) and 100 cm -1 (full triangles). 
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FIG. 4: (a) Small polaron eigenenergies for a nanowire containing AT = 15 molecules. The temperature is fixed to T — 100 K 
and the small polaron binding energy is equal to Eb = 3 cm -1 . The full line represents the density of states of the corresponding 
infinite lattice with translational invariance and characterized by the core hopping constant of the finite size nanowire. (b) 
Eigenfunctions of the two eigenstates with the higher energy. 
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FIG. 5: (a) Small polaron eigenenergies for a nanowire containing N = 15 molecules. The temperature is fixed to T = 100 
K and the small polaron binding energy is equal to Eb = 15 cm -1 . The full line represents the density of states of the 
corresponding infinite lattice with translational invariance and characterized by the core hopping constant of the finite size 
nanowire. (b) Eigenfunctions of the two eigenstates with the higher energy. 
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FIG. 6: (a) Small polaron eigenenergies for a nanowire containing N = 15 molecules. The temperature is fixed to T = 100 
K and the small polaron binding energy is equal to Eb = 7.2 cm - . The full line represents the density of states of the 
corresponding infinite lattice with translational invariance and characterized by the core hopping constant of the finite size 
nanowire. (b) Eigenfunctions of the two eigenstates with the higher energy. 
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FIG. 7: Behavior of the energy of the localized states from the top of the polaron band vs the small polaron binding energy 
for two values of the lattice size = 10 (black symbols) and N = 15 (white symbols). The temperature is fixed to T = 100 K. 
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FIG. 9: Phase digram in the space of the parameters (a, 0) (see the text). The full lines characterize the occurrence of high 
energy localized states located above the polaron band. The lines discriminate between three regimes connected to the presence 
of zero (OHE), one (1HE) or two (2HE) high energy localized states. The dashed lines refer to the occurrence of low energy 
localized states, located below the polaron band, and discriminate between three regimes connected to the presence of zero 
(OLE), one (1LE) or two (2LE) low energy localized states. 
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FIG. 10: Phase digram in the space of the parameters (Eb, T) for TV = 10 (full lines) and TV = 20 (dashed lines). For each N 
value, the diagram exhibits two critical curves which discriminate between phases with zero, one and two high energy localized 
states. 



